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Decaying dark matter with heavy axino 
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A TeV scale decaying dark matter chiral multiplet N is introduced in addition to the minimal 
supersymmetric standard model(MSSM). For a calculable abundance of N, we introduce heavy 
axino decaying to N and MSSM particles including the lightest supersymmetric particle(LSP). In 
the scenario that heavy axino, once dominating the energy density of the universe, decays after 
the LSP decouples, it is possible to estimate the relative cosmic abundances of N and the LSP. 
Dimension 6 interactions allow the lifetime of the fermionic or the bosonic superpartner of N in 
the 10 27 s range to be compatible with the recent astrophysical bounds. A diagrammatic strategy 
obtaining a suppression factor 1/M 2 is also given. 
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I. INTRODUCTION 



Recent observations of high energy galactic positrons, 
electrons, antiprotons, and gamma rays attracted a great 
deal of attention on dark matter(DM) scenarios. If TeV 
scale decaying DM(DDM) decays at the present epoch to 
produce these high energy particles, it needs a very long 
lifetime (of order 10 27 s) so that its decay is within the 
allowed limits of experimental observations. On the the- 
oretical side, the standard modcl(SM) has been extended 
to the MSSM, mainly to solve the gauge hierarchy prob- 
lem. This supersymmetric(SUSY) extension has found 
another bonuses: the existence of the LSP x as a c °ld 
DM candidate and the gauge coupling unification around 
(2 - 3) x 10 16 GeV. In Ref. Q, a further extension of 
the MSSM by an additional DM component N (called 
A^dmMSSM) to produce enough high energy positrons 
was suggested together with charged SM singlets i? ± to 
explain PAMELA'S excess positrons 0] from N + x anni- 
hilation. Interestingly, grand unified theories(GUTs) al- 
lowing charged SM singlets -E are possible in the flipped 
SU(5) GUT [|, which has an ultraviolet completion in 
the heterotic string 0]. 

Later last year, the ATIC data raised the DM scale up 
to TeV @, and the genie for TeV scale DDM has been 
let free. With TeV DDM N, the charged singlet E ± of 
Ref. [lj may or may not be needed below the mass scale 
tojv but TVdmMSSM can still be considered. However, 
the recent Fermi LAT data is in conflict with the ATIC 
data of several hundred GeV electrons @ . Even though 
the TeV scale cosmic-ray(CR) electrons are explained by 
the known astrophysical backgrounds, PAMELA'S CR 
positron excess at the 10-80 GeV range may need an- 
other contribution beyond the known backgrounds hfl. 
On the other hand, PAMELA'S low antiproton flux ||| 
has been generally regarded as a difficulty of DDM sce- 
nario [l], Q ■ Note however that the old background es- 



timates of the H.E.S.S. data [lfj had large systematic 
uncertainties. For example, Ref. [ll[ considered these 
uncertainties to allow a leptonic background smaller by 
a factor 0.85 of the old background value. If one applies 
this argument to antiproton flux also, one can allow some 
antiprotons from DDM decay Interestingly, PAMELA'S 
CR antiproton flux above 10 GeV has the same shape as 
the old background estimate, which may be interpreted 
as "the old estimate in fact contained extra antiprotons" . 
This new explanation of the old antiproton background 
allows a room for antiproton injection to the galaxic DM 
soup from DDM decay. So, models producing some an- 
tiprotons in addition to positrons need not be ruled out 
from the outset. 

With this new perspective, now it is very interesting 
to consider the TeV scale DDM possibility, even allowing 
some antiproton flux from the DDM decay though we will 
skip the discussion of the antiproton flux in this paper. In 
this spirit, we consider the TeV scale DDM possibility by 
the simplest extension of the MSSM with just one chiral 
multiplet N at the next mass level beyond the MSSM, 
which is an AT DM MSSM model [l[. The supermultiplet N 
contains the bosonic partner N and the fermionic partner 
iV. 1 The chiral field N becomes a two-component mas- 
sive Majorana fermion at the true vacuum. The LSP \ 
is assumed to be stable with the unbroken R-parity and 
may constitute a dominant portion of galactic DM. Then, 
the TeV scale DDM TV can decay to MSSM particles. The 
needed range of the N lifetime with the stable LSP \ is 
~ {m-N /m x )l0 26 s. 2 The number density of the N chiral 
multiplet is completely unknown at this point. But, if 



1 Without confusion, we use the same notation N for the super- 
multiplet and its fermionic partner. 

2 The lifetime as a function of mjv/m x for two DM components 
can be gleaned from 
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some heavier particle X dominates the energy density of 
the universe and decays to both N and \ below the LSP 
decoupling temperature, it is possible to estimate the rel- 
ative abundances of N and \- We explore this possibility, 
interpreting X as the axino [l3[. 

The axion has the anomalous coupling to gluons. So, 
the heavy axino enables us to estimate the relative abun- 
dances of N and \ through the anomalous coupling and 
a supcrpotcntial term, 
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where c g are coefficients of 0(1), a g 



are the gauge 
couplings, and $ g s are the vacuum angle terms. $3 
defines the axion: 03^3 = X/F a [1J|. The rele- 
vant axino decay Lagrangian [l5[ is ((X) /M')NNa + 
{az/AV^nFajGa^G^a where G is the gluino and G» v 
is the gluon field strength. Here, we neglect the coupling 
4> u (f>dXX/AM', assuming that the LSP is predominantly 
bino. One gluino will produce one LSP in the end, and 
hence we expect the following N and \ ratio from the 
axino decay, in the limit m^r 2> Tn x , 
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FIG. 1: Possible operators for N decay suppressed with a 
super-heavy mass 1/M 2 : (a) The propagator with one arrow 
accompanies a momentum in the numerators of the propaga- 
tor, and (b) A diagram with the A-term. 
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(2) Majorana mass m n in the numerator of a propagator as 
im n j(k 2 — m„). For Fig. [Ha), we have an operator for 



10" 



To obtain this ratio at the level of ~ 
and F a ~ 4 x 10 11 GeV, we need M' ~ 2 x 10 15 GeV 
which falls in a broad GUT scale with our notation of 
M GUT - 10 15 - 5 x 10 16 GeV. 



II. MODELS 

In addition to the MSSM symmetries we introduce the 
R-parity and the Peccei-Quinn(PQ) symmetry U(l)r- In 
addition, we also introduce matter parity P. The at- 
tractive feature of the PQ symmetry is that it solves the 
strong CP problem, the resulting invisible axion may con- 
stitute a cold DM component, and its breaking scale is 



14] 



narrowed down to a window 3 10 9 < F a < 10 12 GeV 
so that our estimate of the N lifetime is more or less 
predictive. 

The simplest 1/M 2 suppression results with four ex- 
ternal fields which however cannot be expressed as a su- 
pcrpotcntial term. This interaction includes a derivative 
coupling. An example of the derivative interaction with 
four external lines is given in Fig. QJa) with the coupling 
e^EN of Ref. [l|. The Weyl field propagator for one 
direction arrow is ifc E /{k 2 E — m E ) and Fig. [2(a) gives 
a dimension 6 operator with one derivative multiplied 
by l/m E . In contrast, two colliding Weyl fields gives a 



3 But, note that there exists the possibility that F a can be larger 
than 10 12 GeV for a small initial misalignment angle lid . 



N decay, i.e. the interaction Lagrangian becomes 
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where the super-heavy mass is M 2 = m E /fjf*j. By 
the interaction {3}, N (or N) decays to (or N) if N 
(N) is heavier than A^ (A^). The dimension 6 interac- 
tion can also arise with five external lines expressible 
as (f^/M 2 )^ 

a,boson0/3,boson^7,fermioii-^boson-^fermion? 

where M is at a GUT scale, and 4> a ,<t>p and § 1 are 
the MSSM chiral fields. For this interaction, a quintic 
superpotential can be written as ~ -^p-A^JUMSSM where 
W^mssm is the dimension 3 superpotential. 

A simpler form compared to ([3]) can be taken as 



A 



N*e c jNej + h.c. 



(4) 



If A^ of Eq. Q develops a vacuum expectation 
valuc(VEV) V, then we obtain the form (T5|) with A = mj 
and (N*) = V* . In fact, there exists a diagram, lead- 
ing to Eq. ([4]) with the ^4-term insertion, as shown in 
Fig. []Jb). The importance of this diagram depends on 
how SUSY is broken, and it is probably more important 
than Fig. [2(a) m the gravity mediation scenario. We 
find two more diagrams comparable to but not exceeding 
Fig. [TJb): one type with m^/A times Fig. [TJb) and the 
other with the insertion of the S-term Bm E E c E . In this 
paper, we study the form (g]) for a concrete discussion. 

Let us introduce X\ t 2 to break the PQ symmetry [ij} 
by their VEVs. Axino is predominantly X, a combina- 
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TABLE I: Color singlet chiral fields and their quantum num- 
bers. 



tion of X x and X 2 . The U(l) r -SU(3) co i or -SU(3) co i or 
anomaly(rCC anomaly) is present either by the heavy 
quark(s) and/or the SM quark(s). The quantum num- 
bers (the R-parity, the hypcrchargc Y, the PQ charge T, 
and the matter parity P) of the needed fields are listed 
in Table [H including the Higgs doublet pair <f> u and 4>d- 
The coupling f'NnXi gives the first term of Eq. fTJ) with 
M' = m n /f' 2 . The lightest among E a is called E which 
has a huge mass splitting from the other E a s, but here 
E is still treated as super-heavy toward a DDM scenario, 
unlike in Ref. [![■ Ts could be assigned such that fam- 
ilies are not distinguishable, nevertheless in Table [T] we 
devised a scheme toward a much smaller electron mass 
compared to the muon and tau masses. 

Toward the fj, term generation a la Giudice and Masiero 
[lH , we also listed £ in Table |T] for the Kahler potential 
term H*<j) u <j)d/Mp. The following renormalizable super- 
potential is consistent with the symmetries of Table (H for 
I = e, fi, t and a = 1, 2, 3, 

1 



W = fi a N ei E c a + m a E E a E c a + -m, 



1 AXfS (5) 

where the first term is the coupling considered in [l[. 



Lifetime of N — The TeV scale sector contains fields 
N,N, saxion s and axino a, whose mass hierarchy is as- 
sumed to be 



m x < m N ,m^ < m s ,m 3 / 2 ,m ii , 



(0) 



where TO3/2 is the gravitino mass, and we use the follow- 
ing VEV of N and the suppression mass, 



(N) ~ V, Mfj 



fifj 



(7) 



We distinguish the three cases for the interaction ^ ac- 
cording to the masses of the N supermultiplet and the 
VEV of N. The discussion on how N develops a VEV 
or not is outside the scope of this paper, and below we 
will simply choose the cases of V = or V 7^ toward 
a phcnomenological study. In Fig. [2] we show some pa- 
rameters fitted to the observed flux. 

Case (a): mjy > m^,V = 0: The matter parity P is 

unbroken and N is a stable particle. If gravitino is lighter 
than the N multiplet, N can decay to gravitino with the 
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FIG. 2: A few fits of DDM masses to CR e ± with ft = f CE . 
The slepton masses are those given at the benchmark points 
of SPSla [Is|. The skyblue (goldenrod) line is the back- 
ground for Case (a) with 70 = 2.60(2.54) accounting the 
galactic propagation through the GALPROP [20T |. The dan- 
delion line is for 70 = 2.42. Propagation of excess from the 
DDM decay is calculated by the GALPROP code. The solar 
modulation potential relevant near the GeV region is $ = 500 
MV. The data are from the CAPRICE(peach dots) [5l[, 
the AMS(navyblue squares) [2^], the HEAT(maroon squares) 
22 |. the PAMELAe* ratio(midnightblue dots) in the inset 
% the calculated PAMELAe* with 70 = 2.54, 2.60(purple, 
midnightblue circles), H.E.S.S.(blue squares) [10(. PPB- 
BETS(green squares) J3, ATIC(gray dots) M, and the 
Fermi LAT(red dots) 



amplitude suppressed by one power of the Planck mass. 
So, we forbid this by the hierarchy, ([6]). By interaction 
flU) , N decays to N* + ej + e c f , and to their charge con- 
jugated states, 
as 



So, we estimate the decay width of N 
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(l+»-2f) 

= |(1 -x + y) 



(8) 



-in 



where x = m 2 ^/rn 2 N ,y = 
A which is of order m 3 / 2 
y — > 0, we obtain T = . 
((HI) to give an order of 10 26 (m_/v/m x ) 
GeV and A = 1 TeV, we need a relation 
10 15 (mAr/TeV) GeV. 



\/rn 2 N and we take a real 



In the limit of 



/768tt 3 M 4 



/V 



and 
For Eq. 



s for m x = 100 



M 



2.8 x 



Case (b) > mjv, V = 0: The matter parity P is un- 
broken and N is a stable particle. The N decay to grav- 
itino is forbidden by the hierarchy, ((6]). By interaction 
dU) , N decays to N + e j + ef , and to their charge con- 
jugated states. So, we estimate the decay width of N 
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MS 

1 — 2 e 7T a Af 4 J VI ' (1 + 5-277) 2 ' ^ 

?7max = 5(1 + i - y) 

where i = mj^/rn 2 ^ and y = m^/m 2 ^. In the limit of 
£ = y = 0, we have T = (A 2 m 3 / 1536tt 3 M 4 ). For 
Eq. (J9j) to give an order of 2 x 10 26 (m^/m x ) s (an- 
other factor 2 for both N and TV* decays), we need 
M ~ 2.8 x 10 15 (tojv/TcV) GeV for x -> 0,y -> and 
m x = 100 GeV and A = 1 TcV. 

Case (c) TOjy > m^r, V 7^ 0: This is the simplest case. 
The matter parity P is broken by the VEV V. Since 
we introduced only one global symmetry U(l)r, the EW 
scale VEV of the PQ charge carrying field N does not 
lead to any other Goldstone boson in addition to the 
one already introduced at the scale F a . The lightest P 
odd particle N decays to the MSSM particles. DDM is 
the fermion V which can decay by the interaction {4}: 
N c (= N) — > e c j + if, and to their charge conjugated 
states. So, we estimate the decay width as 

which becomes 3xl0- 23 s~ 1 (10 15 GcV/A/) 4 (V/100 GeV) 2 
•(^4/10TeV) 2 (mAr/TeV), where we neglected the slepton 
mass. To give an order of l0 26 (rriN /m x ) s, we need 
a relation M ~ 7.4 x 10 15 GcV(mjyV/ TcV 2 ) 1 / 2 for 
m x = 100 GeV. 



III. FITTING TO CR ELECTRONS AND 
POSITRONS 



curve as shown in the inset 0. The PAMELA data does 
not give an independent flux for (e + + e _ ), and hence we 
calculate the total flux from the ratio, r = e + /(e + + e~), 
using the calculated background estimates of the gold- 
enrod and skyblue curves for Case (a) (2(j. The mid- 
nightblue dots in the inset go to the purple (70 = 2.54) 
and midnightbluc(7o = 2.60) circles for e + + e~ . The 
production rates of e ± from the DDM decay are calcu- 
lated using the isothermal profile. PYTHIA is used to 
obtain the e ± spectrum from the decay of DDM. The 
galactic propagation of these CR (from the DDM de- 
cay and the local sources) to Earth is estimated using 
the CR propagation package GALPROP (2(|, partially 
modifying it. The parameters of the fitted curves are as 
shown in the figure. For example, Case (a) with the ma- 
genta dashed line is an excellent fit with = 6 TeV, 
rrifj = 1 TeV at a benchmark point of the SPSla [l9[ 
with the couplings f e E = and / m e = |/ t b. From the 
figure, we notice that the fit is a combination of 70, the 
DDM mass and the couplings. 



IV. CONCLUSION 

We introduced just one more chiral multiplet N be- 
yond the MSSM particles at the next higher mass level 
of TeV, which allows the N lifetime in the 10 27 s range by 
dimension 6 operators. We have successfully fitted both 
the PAMELA and Fermi LAT data with e ± produced 
by the decay of N . The heavy axino, decaying to both 
N and \ below the x decoupling temperature, enables 
us to estimate the relative abundances of N and x- An 
interesting aspect of this axino decay scenario is that the 
suppression mass scales considered in Eqs. (|1I8I9I10P fall 
in the general GUT scale Mqut- 



In Fig. [21 we present the best fit DDM masses at the 
SPSla benchmark point [lj| for Cases (a, c). Cases (a) 
and (6) with the exchange N <-> TV are almost indistin- 
guishable. The Fermi LAT data may be fitted by a dif- 
ferent injection spectrum 70 = 2.42 (the dandelion line), 
but then the PAMELA data is far above this dandelion 
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